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Thermodynamic functions for ammonia as an ideal gas at one atmosphere pressure have been 
evaluated. The contribution of the highly anharmonic out-of-plane vibrational mode, including its 
large coupling with rotation and its coupling with the other vibrational modes, is considered in detail. 
Tables of C%IR, (H° — Eq)IRT, {E°q — G)IRT, and S°/R have been calculated at closely spaced inter- 
vals from 50 to 5000 °K within an overall uncertainty of less than 0.1 percent at 1000 °K. 
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1. Introduction 

This paper is part of a program originating with the 
Office of Standard Reference Data of the National 
Bureau of Standards to obtain accurate thermody- 
namic properties for ammonia. Subsequent papers of 
this program will include real-gas corrections to the 
ideal-gas properties contained herein. 

In this paper the ideal-gas thermodynamic properties 
of ammonia in the ground electronic state are calcu- 
lated at closely spaced temperature intervals from 
50 to 5000 °K. The calculations include the contribu- 
tions of ordinary vibrational anharmonicity and vi- 
brational-rotational coupling together with rotational 
stretching and rotational quantum effects. Particular 
attention is given to the treatment of the rotational 
and vibrational anharmonic effects due to the molecu- 
lar inversion. To this end a method is developed for 
calculating this contribution such that an overall 
accuracy for the thermodynamic properties of am- 
monia approaching that of the simplest diatomic 
molecules is obtained. 

In the vibrational ground state ammonia is a pyram- 
idal molecule with the nitrogen atom approximately 
0.38 A above the plane of the three hydrogen atoms. 
The distance between the nitrogen atom and each of 
the hydrogen atoms is approximately 1.0 A. Thus the 
pyramid is relatively flat and the potential barrier to 
vibration of the nitrogen atom through the plane of 
the hydrogen atoms so as to invert the pyramid is 
low (approximately 2000 cm -1 [l)). 1 

Ammonia has six vibrational degrees of freedom. 
However, owing to symmetry there are only four normal 
vibrational modes, two of which have double degen- 
eracies. The vibrational structure may thus be char- 
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acterized by the six quantum numbers (vitfefja&i*)) 
where the v\ are the principal vibrational quantum 
numbers associated with the vibrational modes v\ 
and the / characterize the small coupling of the de- 
generate i>3 and Va asymmetric vibrational modes with 
the angular momentum of the molecule. The V\ and 
v>> modes are symmetric, nondegenerate vibrations 
corresponding to stretching the N-H bond and spread- 
ing the pyramid, respectively. (The v* mode is also 
referred to as the out-of-plane vibration.) The rotational 
structure for the ground state is that of a symmetric 
top, so that the rotational energy levels are charac- 
terized by the two quantum numbers ,/ and K. 

The v-z mode is the vibration associated with the 
inversion of the pyramid. Each vibrational state of this 
mode is split into two components, corresponding 
respectively to symmetric and antisymmetric eigen- 
functions. (In the nomenclature these are identified 
by the superscript "s" or "a" to the quantum number 
V2>) The spacing of these levels is highly anharmonic 
except for vibrational energies sufficiently higher or 
lower than the potential barrier to inversion. Thus the 
vibrational levels for the ground state of the V2 mode 
are only slightly split, of the order of a fraction of a 
wavenumber, but the splitting increases rapidly with 
increasing vibrational quantum number. At energy 
levels sufficiently elevated above the barrier (greater 
than about 3000 cm -1 ) the levels become uniformly 
spaced with a spacing of approximately half that of 
the V2 fundamental. Similarly, due to inversion the 
rotational structure undergoes a transition from that 
of a pyramid for the low values of the V2 quantum num- 
ber to that of a planar molecule for the higher values, 
so that the rotational motion is also strongly coupled 
to the V2 mode. For a more detailed discussion of these 
effects, see the review by Herzberg [2| of the pyramidal 
XY S molecule and Nielsen [3|. The energy structure 
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of the other vibrational modes is similar to that of 
typical, slightly anharmonic rotating vibrators. 

In section 2 the equations for the calculations are 
developed. This is followed in section 3 by a discussion 
of the molecular data, and in section 4 by a discussion 
of the calculated properties and an evaluation of their 
accuracy. Finally, a discussion is given in section 5 
of some of the earlier work. 

2. Partition Function 

It is convenient to represent the rotational and vi- 
brational structure by an expression for which the 
energy levels of the v 2 mode and its coupling with 
the other degrees of freedom are written explicitly. 
Thus the energy levels in wavenumbers are written 

T(v,J 9 K)=G'(v)+F'(v,J 9 K) + Gi(v 9 J 9 K), (1) 

where the primes denote that dependence on v 2 
is omitted. The first term in eq (1) describes the 
vibrational energy characterized by (v\0v l ^v l 4 4 ) and in 
notation similar to that in Herzberg [2] is approximated 
by 

G'(v)—^ i vm+'^XijViVj + '^ i xu(vi — l)vi +. . . 



where the v% are fundamentals and the x and g are 
relatively small anharmonicities. The second term in 
eq (1) is approximated by 

F' (v, 7, K) = B'JU + 1) + {C'-B' V W (3) 

where C v is the rotational constant corresponding 
to rotation about the symmetry axis and B v is the con- 
stant for axes perpendicular to the symmetry axis. 
The quantities D J V , D J V K , and D§ are small coefficients 
associated with centrifugal distortion of the mole- 
cule. The subscript v to each of the coefficients of 
eq (3) indicate dependence or coupling with the 
vibrational motion. The primes indicate that de- 
pendence on vi is omitted; that is, the rotational 
constants of eq (3) are evaluated for v 2 equal to zero. 
The coupling of the rotational motion to vibrations 
other than the Vi mode is relatively small and is 
included here only to first order in the fundamentals. 

B' v = B -^'afvi + 

i 

C' v = Co — 2' cttyi -h , i 7^2 

i 

where B and C are rotational constants, and the a 
are small corrections. All coordinates involving 
the v 2 mode, including its coupling with the other 



degrees of freedom are included in G 2 . Thus G 2 
includes those effects characterized by 

{v x v 2 v!?v{*JK) - {vfivifv^JK). 

The partition function for the rotational and vi- 
brational structure may be written 



(?=2 exp{-T(v,J,K)c 2 IT}, 

v,J,K 



where T(v, 7, K) is discussed in eqs (1), (2), and 
(3). c 2 is the second radiation constant and T is the 
absolute temperature. The sum extends over all 
rotational and vibrational degrees of freedom in- 
cluding all values of the v, J, K quantum numbers. 
In the evaluation of eq (4) the contributions of the 
vibrational anharmonicity and vibration-rotation 
coupling effects contained in eqs (2) and (3) are 
all to be accounted for to first order; that is, only 
terms linear in x, g, and a are retained. Though 
G 2 is not conveniently described by the usual power 
series in the vibrational quantum numbers, its con- 
tribution, including coupling effects with the other 
degrees of freedom, is also to be included to this order. 

The contribution of the G 2 term is obtained by ex- 
panding the partition function about the term-by- 
term sum of the observed (Oi^OO) energy levels. This 
summation incorporates the major effect of the highly 
anharmonic v 2 mode, which is the coupling between 
the symmetric and antisymmetric states, but it does 
not include the contribution of the v 2 coupling with the 
other vibrational modes. From an examination of the 
spectroscopic vibrational data (see Benedict et al. 
[4]), it can be seen that this cross coupling produces a 
shifting of energy levels of about the same magnitude 
as that of the other vibrational couplings (those repre- 
sented by the x and g terms in eq (2)). Thus the con- 
tribution to the partition function of this cross coupling 
should be of the same order as that from the other 
anharmonicities, and a first order approximation for 
this contribution should be adequate. 

In the expansion that follows expressions for the 
cross coupling are developed and their contribution 
for each individual (Oi^OO) state obtained explicitly. 
The cross coupling between a pair of vibrational states, 
e.g., states corresponding to the particular quantum 
numbers v\ and v 2 , is defined as 



i G(viV 2 00) - G(v,000) - G(0v 2 00) , 



(5) 



(4) 



G Vll)2 = k l2 G (vtOOO) -G(0v 2 00) + . 



where G{v\V 2 v§v l f) is the vibrational energy. A useful y 
representation for the coupling involving the highly 
anharmonic v 2 mode is the development of G VlV2 as a 
power series of products of the observed energies, 
that is 



(6) 



The constant k\ 2 is readily determined in terms of the 
observed coupling. Thus for coupling between the 
fundamentals the observed coupling x\ 2 is given by 
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*i2 = A;i2G(1000) -r;(0100) + . . . 

Note that eq (6) reduces identically to the power series 
of eq (2) when the cross coupling is between nearly 
harmonic modes. 

Consider the energy of the nearly harmonic oscil- 
lator v m , neglecting all anharmonicities but the cross 
coupling with v 2 . The energy levels for such an oscil- 
lator are given by 

Substituting from eq (6) and defining E 2 = G(Qv 2 Q0), 
tli is becomes 



The prime indicates that the m = 2 factor in the product 
corresponding to self coupling of the v 2 mode is 
omitted. 

The coupling of rotation with the v 2 mode is in- 
cluded by evaluating the classical rotational partition 
function for each of the ((feOO) energy levels and em- 
ploying these as separate weighting factors to each of 
the terms in eq (12). These factors are the usual ex- 
pression for a classical nonrigid rotator, 

(?/ i (^)=7r 1 / 2 (a-|crc)- 1 ^ 2 (l + p (1) r-r-p( 2 )P). (13) 

The quantities cr# and 07 are given by, 



(Th z 






G n 



V m Vm + k-imV 



•ImV mi'm 



E,_. 



(7) 



The contribution of G m to the partition function is the 
factor 



(8) 



Qm = \ ^ ex P ~~ [( c 2/r) {VrnVm + fc 2 m#2*'ml>m)] 



where d m is the degeneracy of the rath mode. The 
anharmonicity term in the exponential of eq (8) is now 
expanded and only the linear term retained to yield 



Using the identity 



T 



V m Vm I J 1 —j;k2mV m V m E 2 \ 



(9) 



- — V 

da^ 



ve 



the contribution of Q m becomes 

-dm 



Q» 



1 — e t 



I-/2, 



(10) 



where 



0"C rp v>r 2 -» 



where B v . z and C V2 are the observed rotational con- 
stants for the particular (Oi^OO) energy level. The factor 
1 -\- p {1) T + p (2) T 2 is the contribution of rotational 
stretching (see, for example (5 and 14)), where the p 
are calculated from the rotational stretching coeffi- 
cients that apply to the particular (Oi^OO) level. 
One further simplification is required in order to 
evaluate the out-of-plane contribution. This is to 
obtain an approximation for the contribution of the 
energy levels corresponding to high values of the v 2 
quantum number. Note that the v 2 mode is nearly 
harmonic for energies corresponding to the fourth 
quantum number and higher, so that the contribution 
of these levels may be replaced by the well-known 
expression for a slightly anharmonic rotating vibrator. 
Combining eqs (12) and (13) and including the closed 
form approximation for the contribution of the states 
corresponding to quantum numbers given by v 2 > 3 
yields the expression for 62 (v, ./, K) used in this 
calculation. 



Bo /CoV /2 



V 2 <4 



&* = j2^lgj d+P (i) ^+P (2) n. 



IT (W«)e~ 



c-i 



hm = d m Y hmEiVml (e c ' 2 



civmlT _ 



(ID 



The factor (1— /2m) is the contribution of the inter- 
action of a particular (6V2OO) energy state with the 
mth vibrational mode, and it depends explicitly on 
Ei. It is clear then that the effect of this coupling is 
to modify each of the terms in the term-by-term sum 
over the (Ot^OO) energy levels by the factor (1 —fimj* 
The term-by-term sum including cross coupling with 
the other vibrational modes is thus given by 

Q v% (vib) = 2 n ( J -/**) ex P (-<*B*IT). (12) 



1/2 



exp (-c 2 E 4 JT) B^/Co 
1-exp (-CtvilT) ' B 4 \CJ 

na-^Mi + cMO + p^), 



(14) 



where E = E 2 and the subscript to E now refers to the 
particular v 2 quantum number. The quantity in the 
braces in eq (14) includes levels only up to v 2 — 3 a . 
Note that the rigid rotational contribution for the 
ground state has been factored from eq (14). In the 

exp (—c 2 E 4 s/T) 



second term the factor 



is the har- 



1 — exp (—C2P2IT) 
monic oscillator contribution with E 4 * as the first level 
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and v'y the fundamental. The factor 1 — 2 . ; is the ordi- 
nary anharmonicity coupling with the y'th vibrational 
mode, where 



cz 



Xij(lj(tj 



iexp[c i v,IT)-l)(exp(c i v J IT)-iy 



e„- 



[ (xu-giil3)di(di+ 1) 
(exp (c 2 ^,/r)-l) 2 



(15a) 



The factor 1 + (j>> is the well-known contribution of 
rotation-vibration coupling, where 



( ofi 1 oiF 



(15b) 



Details concerning the derivation of quantities in eqs 
(15) may be found in Friedman and Haar [6], Woolley 
[7], and Pennington and Kobe [8J, for example. 

The contributions to the partition function of the 
two remaining terms G' (t>, J, K) and F'(v, J, K) are 
quite straightforward (see [6, 7, 8]). With these the 
complete partition function is finally obtained. 



1SJ 

n V 2 



.°"«..°'f'o 



Q* ■ <2tr " Q V . (16) 



As noted earlier the primes on the product summations 
indicate that factors containing the vibrational coordi- 
nate vi are omitted. The quantity S is the symmetry 
number. For the usual symmetric top molecule this has 
the value three. However, each of the components of 
the nondegenerate but split (Oi^OO) levels have been 
counted with a statistical weight of unity, so that the 
factor Q V2 is larger by a factor of two. This error is cor- 
rected by setting S = 6. (Incidentally S = 6 corre- 
sponds to the symmetry for a planar configuration for 
NH 3 , and eq (16) might be thought of as an expansion 
about this planar configuration.) The product over m is 
the harmonic oscillator contribution for the modes 
other than v<l\ the products over ij and over k are the 
corresponding anharmonicity and rotation-vibration 
coupling contributions, respectively. The degeneracy 
coefficients d are unity for m= 1, 2, and two for m=3, 4. 
The Oij and </>a- are given by eq (15). Q tT is the usual 
contribution of translation. Lastly, the factor Q qr is a 
correction for rotational quantum effects. It is very 
small at temperatures for which the vibrations are 
appreciably excited, so that it is not coupled to the 
i>2 mode. The contribution Q qr is given by Stripp and 
Kirkwood [9] for an asymmetric rotator, which for a 
symmetric top rotator reduces to 

^ 1+ 3\ ] 4o-J 15 V 1 4o-c 32o-*J 



3. Molecular Data 

The molecular data upon which the calculations are 
based are listed in tables 1 and 2. Table 1 lists the data 
sensitive to the inversion coordinates and utilized in 
the term-by-term sum (the quantity in braces in eq (14)). 
The vibrational energies corresponding to V2 = 0, 3 are 
from [4]. The energies corresponding to tte=l, 2 are 
from Garing et al. [10]. The values listed for the rota- 
tional constants B v and C v corresponding to the 
V2 — states are from [4]; those corresponding to the 
V2= 1, 2 states from [10]; and to the #2 — 3 states from 
Benedict et al. [11]. The rotational stretching constants 
D J V , Dlf, D$ are the average of the symmetric and anti- 
symmetric values. The values for V2 — are from [4]; 
for i>2=l, from [10]; for #2 — 3, from [11]; finally, data 
for #2 = 2 were taken equal to the V2 = 1 values. 



TABLE 1. The out-of-plane molecular constants in cm -1 , 
v 2 «£ 3, (0v 2 00) 



State 


0* 0" 


I s 1" 


2* 2° 


3« 3" 


Vt 


0.0 0.793 


932.41 968.04 


1597.42 1882.16 


2383.46 2895.48 


Bv 2 


9.9443 9.939 


10.07 9.89 


10.26 9.73 


9.50 9.20 


c Vt 


6.196 6.198 


6.05 6.13 


5.90 6.13 


6.16 6.26 


D i 


8.1 X10- 4 


8.9 xlO- 4 


8.9 X 10- 4 


-3.3 xlO- 4 


QJK 


- 14.5 x 10- 4 


-17.5X10- 4 


-17.5 XlO- 4 


15.0 X10- 4 


£>* 


7.5 X10- 4 


10.5 XlO- 4 


10.5 x 10- 4 


-11.3 XlO" 4 



Table 2. 



Vibrational, rotational and coupling constants in cm- 
see eqs (2), (3) 



p,(1 0) 


p! 2 = E 4a — E 4s 


M0 o l 0) 


v 4 (0 1) 


3323.06 


475 


3443.381 


1627.77 


xu=-30 


Jt22 = + 9 


x 3 3 = -18.5 


x 44 = -8.8 


x,2 = +20.6 


jt 23 = + 32.36 


*34 = -17.3 


^44 = + 2.65 


* 13 = -92 


JC 2 4 = -10.7 






x, 4 = -6.7 


#4 = 9.35, C 4 = 6.21 






of = 0.135 


of =+0.015 


of = + 0.176 


a? = -0.230 


af= .078 


a 2 c = + 0.098 


a? = -0.009 


a$ = + 0.066 



The remaining molecular data are given in table 2. 
The vibrational fundamental v\ is from Benedict 
et al. [12]; v 3 is from [11]; v 4 is from [10]. The constants 
z^and #22 which correspond to the high energy, nearly- 
harmonic region of the v>i mode are derived from the 
reported values for the (04*00) and (05 s 00) levels at 
3330 cm -1 and 4315 cm -1 given by Walsh and Warsop 
[13] and from theoretical considerations that at high 
energies the spacing between the levels should become 
nearly equal to one-half the (0100) value. Also pertain- 
ing to this region are the rotational constants B 4 and C 4 
which are taken equal to the averages of the values for 
f 2 = 3 s and 3 a . The anharmonicity coefficients #12, #13, 
£i4, #23, #24, #33, and #34 are from [4], and the coeffi- 
cients #44 and g44 are from [10]. The anharmonicity 
coefficient xu is estimated from data for similar 
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stretching vibrations. The coupling constants af 
are from [4|. The small coefficient £33 is taken as zero. 
The values used for the atomic masses are from 
Cameron and Wichers [35J. The values used for the 
other physical constants are given in [36]. 

4. Thermodynamic Functions 

The thermodynamic properties of ammonia as an 
ideal gas at a pressure of one atmosphere are listed in 
► table 3. These include the Gibbs (free) energy function, 
enthalpy function, heat capacity at constant pressure, 
and entropy tabulated in dimensionless units at closely 
spaced temperature intervals from 50 to 5000 °K. 
Should it be desired to obtain the properties for a pres- 
sure of 1 bar (10 6 dynes cm -2 or 10 5 Nm~ 2 ) the magni- 
tude of the tabulated free energy and entropy values 
should be increased by 0.0132. Conversion factors to 
other units in frequent use are given in table 4. The 
entries in table 3 are truncated at four decimals (ex- 



cept for C° p /R above 2500 °K which is truncated at 
three) but at best they are accurate to three decimals. 
It is felt that the additional significance could be useful 
in applications involving differences between the 
entries. The density of entries permits use of linear 
interpolation throughout the table with a resulting 
interpolation error smaller than the uncertainty in the 
tabulated values. 

As indicated by the earlier discussion the present 
calculation, except for the out-of-plane vibrations, is 
based on an expansion about the harmonic oscillator- 
rigid rotator model. In this expansion only first order 
(linear) terms in the coupling between the various 
degrees of freedom are retained. In addition to the 
errors inherent in the derivation of the partition func- 
tion, eq (16), there are also errors due to uncertainties 
in the molecular constants, errors which arise from 
certain resonances between several of the low-lying 
vibrational states, the effects of which are omitted, and 
lastly errors arising from omission of excited electronic 
states. 






Table 3. Thermodynamic functions for ammonia 



The tabulated values are spaced sufficiently close to permit linear interpolation between the entries with a resulting error less than 
their overall uncertainty. For convenience the values at 273.15 and 298.15 °K have also been included. 



T 


(E° -G)IRT 


(H°-E$IRT 


C° P IR 


S°/R 




50. 


12.0017 


3.9539 


4.0029 


15.9555 




60. 


12.7234 


3.9620 


4.0026 


16.6854 




70. 


13.3346 


3.9678 


4.0026 


17.3024 




80. 


13.8648 


3.9721 


4.0027 


17.8369 




90. 


14.3328 


3.9756 


4.0029 


18.3084 




100. 


14.7519 


3.9783 


4.0033 


18.7302 




110. 


15.1311 


3.9806 


4.0039 


19.1117 




120. 


15.4776 


3.9826 


4.0051 


19.4602 




130. 


15.7964 


3.9844 


4.0069 


19.7808 




140. 


16.0918 


3.9861 


4.0097 


20.0779 




150. 


16.3668 


3.9878 


4.0138 


20.3546 




160. 


16.6243 


3.9896 


4.0194 


20.6138 


• 


170. 


16.8662 


3.9915 


4.0267 


20.8577 




180. 


17.0944 


3.9937 


4.0358 


21.0881 




190. 


17.3104 


3.9962 


4.0468 


21.3066 




200. 


17.5154 


3.9991 


4.0599 


21.5145 




210. 


17.7106 


4.0023 


4.0749 


21.7129 




220. 


17.8969 


4.0060 


4.0920 


21.9029 




230. 


18.0751 


4.0101 


4.1110 


22.0852 




240. 


18.2458 


4.0147 


4.1319 


22.2606 




250. 


18.4098 


4.0199 


4.1548 


22.4297 




260. 


18.5676 


4.0255 


4.1794 


22.5931 




270. 


18.7196 


4.0317 


4.2057 


22.7513 




273.15 


18.7663 


4.0338 


4.2143 


22.8001 




280. 


18.8664 


4.0384 


4.2336 


22.9048 




290. 


19.0082 


4.0456 


4.2630 


23.0539 




298.15 


19.1204 


4.0520 


4.2882 


23.1724 




300. 


19.1455 


4.0534 


4.2939 


23.1989 




310. 


19.2785 


4.0617 


4.3260 


23.3402 




320. 


19.4076 


4.0704 


4.3593 


23.4781 




330. 


19.5330 


4.0797 


4.3937 


23.6127 




340. 


19.6550 


4.0895 


4.4291 


23.7444 




350. 


19.7737 


4.0997 


4.4654 


23.8733 




360. 


19.8893 


4.1104 


4.5025 


23.9997 




370. 


20.0021 


4.1215 


4.5402 


24.1235 
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Table 3. Thermodynamic functions for ammonia — Continued 




T 


(E° Q -G)/RT 


(H°-El)JRT 


c° p ir 


S°IR 


380. 


20.1121 


4.1330 


4.5785 


24.2451 


390. 


20.2196 


4.1449 


4.6173 


24.3646 


400. 


20.3247 


4.1572 


4.6565 


24.4819 


410. 


20.4275 


4.1699 


4.6961 


24.5974 


420. 


20.5282 


4.1829 


4.7359 


24.7111 


430. 


20.6268 


4.1962 


4.7759 


24.8230 


440. 


20.7234 


4.2098 


4.8161 


24.9332 


450. 


20.8181 


4.2237 


4.8563 


25.0419 


460. 


20.9111 


4.2379 


4.8966 


25.1491 


470. 


21.0024 


4.2524 


4.9368 


25.2548 


480. 


21.0921 


4.2671 


4.9771 


25.3592 


490. 


21.1802 


4.2820 


5.0172 


25.4622 


500. 


21.2669 


4.2971 


5.0572 


25.5640 


520. 


21.4360 


4.3278 


5.1368 


25.7639 


540. 


21.6000 


4.3593 


5.2158 


25.9592 


560. 


21.7591 


4.3912 


5.2939 


26.1503 


580. 


21.9137 


4.4237 


5.3712 


26.3374 


600. 


22.0643 


4.4566 


5.4475 


26.5208 


620. 


22.2109 


4.4898 


5.5230 


26.7007 


640. 


22.3540 


4.5232 


5.5975 


26.8772 


660. 


22.4937 


4.5569 


5.6711 


27.0506 


680. 


22.6302 


4.5907 


5.7438 


27.2210 


700. 


22.7638 


4.6247 


5.8157 


27.3885 


720. 


22.8946 


4.6588 


5.8867 


27.5533 


740. 


23.0227 


4.6929 


5.9568 


27.7156 


760. 


23.1483 


4.7271 


6.0261 


27.8754 


780. 


23.2715 


4.7613 


6.0946 


28.0328 


800. 


23.3925 


4.7954 


6.1622 


28.1879 


820. 


23.5113 


4.8296 


6.2291 


28.3409 


840. 


23.6281 


4.8637 


6.2952 


28.4918 


860. 


23.7430 


4.8978 


6.3604 


28.6407 


880. 


23.8559 


4.9317 


6.4249 


28.7877 


900. 


23.9672 


4.9656 


6.4885 


28.9328 


920. 


24.0767 


4.9994 


6.5513 


29.0761 


940. 


24.1845 


5.0331 


6.6133 


29.2176 


960. 


24.2909 


5.0667 


6.6745 


29.3575 


980. 


24.3957 


5.1001 


6.7348 


29.4958 


1000. 


24.4990 


5.1334 


6.7943 


29.6324 


1020. 


24.6010 


5.1665 


6.8530 


29.7676 


1040. 


24.7017 


5.1995 


6.9108 


29.9012 


1060. 


24.8010 


5.2323 


6.9678 


30.0334 


1080. 


24.8991 


5.2650 


7.0239 


30.1641 


1100. 


24.9960 


5.2975 


7.0792 


30.2935 


1120. 


25.0918 


5.3298 


7.1335 


30.4216 


1140. 


25.1864 


5.3619 


7.1870 


30.5483 


1160. 


25.2799 


5.3938 


7.2397 


30.6738 


1180. 


25.3724 


5.4256 


7.2914 


30.7980 


1200. 


25.4639 


5.4571 


7.3424 


30.9209 


1220. 


25.5543 


5.4884 


7.3924 


31.0427 


1240. 


25.6438 


5.5195 


7.4416 


31.1633 


1260. 


25.7324 


5.5504 


7.4900 


31.2828 


1280. 


25.8200 


5.5811 


7.5375 


31.4011 


1300. 


25.9068 


5.6115 


7.5841 


31.5183 


1320. 


25.9927 


5.6418 


7.6299 


31.6345 


1340. 


26.0778 


5.6718 


7.6749 


31.7495 


1360. 


26.1620 


5.7016 


7.7190 


31.8636 


1380. 


26.2455 


5.7311 


7.7624 


31.9766 


1400. 


26.3281 


5.7604 


7.8049 


32.0886 


1420. 


26.4101 


5.7895 


7.8466 


32.1996 


1440. 


26.4912 


5.8184 


7.8875 


32.3096 
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Table 3. Thermodynamic functions for ammonia — Continued 




T 


(El-G)IRT 


(H°-E$fRT 


QIR 


S°IR 




1460. 


26.5717 


5.8470 


7.9277 


32.4187 




1480. 


26.6514 


5.8754 


7.9671 


32.5268 




1500. 


26.7305 


5.9035 


8.0058 


32.6340 




1550. 


26.9252 


5.9729 


8.0992 


32.8981 


> 


1600. 


27.1159 


6.0407 


8.1881 


33.1566 




1650. 


27.3028 


6.1071 


8.2728 


33.4099 




1700. 


27.4861 


6.1720 


8.3532 


33.6581 


> 


1750. 


27.6659 


6.2354 


8.4298 


33.9013 




1800. 


27.8424 


6.2974 


8.5027 


34.1398 




1850. 


28.0158 


6.3579 


8.5719 


34.3737 




1900. 


28.1862 


6.4171 


8.6378 


34.6032 


> 


1950. 


28.3536 


6.4748 


8.7004 


34.8284 




2000. 


28.5182 


6.5312 


8.7599 


35.0494 




2050. 


28.6802 


6.5863 


8.8164 


35.2664 


x 


2100. 


28.8395 


6.6400 


8.8702 


35.4795 


.-- 


2150. 


28.9964 


6.6925 


8.9214 


35.6889 




2200. 


29.1508 


6.7437 


8.9701 


35.8945 




2250. 


29.3030 


6.7937 


9.0162 


36.0966 




2300. 


29.4528 


6.8425 


9.0602 


36.2953 


~, 


2350. 


29.6005 


6.8901 


9.1021 


36.4906 




2400. 


29.7460 


6.9366 


9.1418 


36.6826 


> 


2450. 


29.8895 


6.9820 


9.1797 


36.8715 


2500. 


30.0310 


7.0263 


9.216 


37.0573 




2550. 


30.1706 


7.0696 


9.250 


37.2402 




2600. 


30.3083 


7.1118 


9.282 


37.4201 




2650. 


30.4441 


7.1531 


9.313 


37.5972 




2700. 


30.5782 


7.1933 


9.342 


37.7716 




2750. 


30.7106 


7.2327 


9.370 


37.9432 


^ 


2800. 


30.8413 


7.2711 


9.397 


38.1123 




2850. 


30.9703 


7.3086 


9.422 


38.2789 




2900. 


31.0977 


7.3452 


9.445 


38.4429 


> 


2950. 


31.2236 


7.3810 


9.468 


38.6046 




3000. 


31.3479 


7.4160 


9.490 


38.7639 




3100. 


31.5922 


7.4835 


9.529 


39.0757 




3200. 


31.8308 


7.5480 


9.564 


39.3788 




3300. 


32.0640 


7.6096 


9.596 


39.6736 




3400. 


32.2921 


7.6684 


9.623 


39.9605 




3500. 


32.5152 


7.7246 


9.648 


40.2398 




3600. 


32.7336 


7.7783 


9.670 


40.5119 


H 


3700. 


32.9474 


7.8297 


9.688 


40.7771 




3800. 


33.1568 


7.8788 


9.704 


41.0357 




3900. 


33.3621 


7.9258 


9.717 


41.2879 


> 


4000. 


33.5633 


7.9707 


9.728 


41.5341 




4100. 


33.7607 


8.0137 


9.737 


41.7744 




4200. 


33.9543 


8.0548 


9.743 


42.0091 


► 


4300. 


34.1443 


8.0941 


9.748 


42.2384 


4400. 


34.3308 


8.1317 


9.750 


42.4626 




4500. 


34.5140 


8.1677 


9.751 


42.6817 




4600. 


34.6939 


8.2021 


9.749 


42.8960 




4700. 


34.8706 


8.2350 


9.746 


43.1056 


P* 


4800. 


35.0443 


8.2664 


9.741 


43.3107 




4900. 


35.2151 


8.2964 


9.734 


43.5115 




5000. 


35.3833 


8.3253 


9.725 


43.7086 



The expansion about the harmonic oscillator limit i n ., , v . -> , r , 

.. ii i . . i ,i mode are small provided ™ is large, and are of the 

is essentially a low temperature expansion, and the ^ T fe 

r errors resulting therefrom for a particular vibrational order of about the square of the linear term. The 
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Table 4 



Conversion factors 

to convert tabulated value 

to quantity having dimensions 

indicated below 


Multiply by 


Cal mole- 1 °K- 1 
Cal g" 1 K-! 
Joules mole -1 °K _1 
Joules g" 1 K-» 


1.98717 
0.116682 
8.3143 
0.488197 



largest error of this type occurs for the v± mode. Thus 
for the most sensitive function, the heat capacity, the 
uncertainty at 2000 °K is less than 10~ 2 dimensionless 
units, but increases to the order of 0.1 at 5000 °K. The 
higher order anharmonicities would contribute an un- 
certainty perhaps as large. Though the anharmonicity 
coefficients X\\ and X\% have large uncertainties, the 
errors introduced by these are mitigated by the fact 
that the fundamentals v\ and *> 3 are large. A 50 percent 
uncertainty in %n and X13 results in an error of about 
0.05 dimensionless units in Cp/R at 5000 °K. Setting 
#33 equal to zero introduces negligible error since this 
coefficient introduces a small correction to an already 
small anharmonic contribution (owing to the large 
value of the v^, fundamental). The remaining anhar- 
monicities (other than the #22 value in table 2) are 
probably accurate to better than several percent, and 
the uncertainties resulting from these are compara- 
tively small. Uncertainties arising from errors in the 
rotational fundamentals are also comparatively small, 
as are uncertainties resulting from the use of estimates 
for some of the rotational stretching constants. The 
uncertainties that arise from errors in the rotation- 
vibration coupling constants are perhaps half as large 
as those arising from errors in the anharmonicity 
coefficients. 

The uncertainties in the values used for the funda- 
mental frequencies are less than several tenths of 
a wavenumber, and the uncertainty from this source 
is negligible. Also negligible are the uncertainties due 
to errors in the observed overtones of the v 2 funda- 
mental. 

It is somewhat more difficult to determine the un- 
certainty resulting from estimates used for the mo- 
lecular constants in the second term in eq (14). The 
largest uncertainty here arises from possible errors in 
v' 2 and £4 . An uncertainty of several percent in these 
results in an uncertainty in Cp/R of less than 0.002 
dimensionless units at 1000 °K, which increases to a 
maximum near 2000 °K of less than 0.01 dimensionless 
units. The errors from estimates of the vibrational 
coupling constants used in this term are completely 
negligible below 2000 °K. Even if a 100 percent uncer- 
tainty in their contribution is assumed, the resulting 
uncertainty at 5000 °K still would be less than 0.1 
dimensionless units in C°JR. 

The spectrum of ammonia is further complicated 
by several large resonances that occur between some 
of the low lying vibrational bands, such as the reso- 
nance between the combination bands (2v 4 4- Ps) and 



(v\ + V3) and between the bands 2v 4 and v\ as discussed 
in [4]. Since the resonances tend to displace the re- 
spective bands symmetrically about their unperturbed 
values, it can be shown that the errors are less than 
the neglected higher order anharmonicities. A similar 
type of error occurs in the treatment of the coupling 
of the V2 mode with the other modes. The coupling con- 
stants as reported are the average of the symmetric and 
antisymmetric values. Thus the effect of this approxi- 
mation is to introduce errors similar to those arising 
from neglect of the resonances discussed above. The 
treatment of centrifugal distortion, which includes a 
quadratic term in the temperature, is consistent with 
the other approximations and the error from this 
source at 5000 °K is probably less than 0.05 in Cp/R. 
At the low temperature end there is a small error 
arising from the rotational quantum approximation Q qr 
(eq (17)). At 100 °K this error is less than 10" 3 in C° P /R, 
but is appreciably larger at 50 °K. 

Finally, the U.V. absorption spectrum indicates 
bands extending from ~ 46,000 cm -1 to over 100,000 
cm -1 , as discussed, for example, in [13]. The error 
from omission of these even at the very highest tem- 
peratures considered here is completely negligible. 
It is felt, therefore, that the uncertainty in the tables 
in the temperature region 100 to 1000 °K is probably 
less than 0.003 dimensionless units in Cp/R. This 
uncertainty increases to 0.02 at 2000 °K and then in- 
creases very rapidly to perhaps 0.3 at 5000 °K. A more 
complete list of the uncertainty estimates, including 
those for the enthalpy and entropy, are given in table 5. 

Table 5. Estimates of uncertainty in dimensionless units 



T 


&C° P IR 


A(H°-E° )IRT 


AS°/R 


100 


0.003 


0.002 


0.003 


1000 


.003 


.002 


.003 


2000 


.02 


.01 


.02 


3000 


.05 


.02 


.03 


4000 


.1 


.05 


.07 


5000 


.3 


.1 


.2 



5. Discussion 

Ideal-gas calculations for NH 3 have been reported 
in the literature since the early thirties, but many oft 
these calculations were based on relatively incomplete 
and inaccurate molecular data and on simplified struc- 
tural models. Also, for the most part, the temperature 
ranges were somewhat limited. The references v 
[17-27] are representative of these early efforts. 

Some of the more recent calculations which have 
been widely used include those of Harrison and Kobe* 
[28] who calculated the heat capacity, enthalpy, and 
entropy from 273.16 to 1500 °K. They included the 
major effects associated with the molecular inversion 
by summing the lower states of the V2 mode term by < 
term, but omitted all vibration and rotational coupling 
effects other than rotational stretching. The review by 
Davies reported in the monograph edited by Din [29] / 



214 



includes the specific heat, entropy, and enthalpy from 
200 to 1000 °K. The calculation follows essentially the 
procedure in [28J, except that centrifugal stretching 

> is omitted. Also it appears that the contribution of the 
states in the term-by-term sum of the v>i mode was 
omitted above the level 3 a . The most detailed of the 
previous calculations is that reported by Yungman 

r et al. [30|, (which also appears in [16|), in which the free 
energy, entropy, and enthalpy are tabulated at 100 °K 
intervals from 298.15 to 6000 °K. In addition to the 
term-by-term sum for the vi mode, this work also in- 
cludes several of the first order vibrational and rota- 
tional coupling contributions. It omits, however, the 
contributions due to coupling of the V2 mode with the 
other vibrational modes, and includes a somewhat 
simplified coupling of the vi mode with rotation. These 
calculations are based on essentially the same mo- 
lecular data as those reported here, except that the 
values used for the vibrational coupling constants xw 
and X44 are about twice as large. One of the most widely 
used set of tables is that of JANAF [31]. The reference 
refers to the latest revision for ammonia, September 30, 
1965. This calculation follows the procedures of [28] 
and [30J, though the specific details are not given in 
the text. Finally, there are several recent calculations 
based on somewhat simplified molecular models, 
which are included in references [32-34]. 

The figures 1 and 2 include a comparison of some 
of the previous calculations with those made here. The 
ordinates are the deviations from the present calcula- 
tions, that is the present calculation minus the others. 
Figure 1 is a plot of entropy differences, and figure 2 
is a plot of heat capacity differences. It is most inter- 
esting to compare the data of Yungman et al. [30], 




FIGURE 1. Dimensionless entropy differences. 
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curve #5 in figure 1. At the low temperatures the 
agreement is quite good, but as the temperature is 
increased, a positive deviation occurs until about 
3000 °K, when the deviation becomes strongly negative. 
The difference at the lower temperatures is probably 
due to their simplified treatment of the coupling of the 
v-y mode with rotation. At temperatures above 1000 °K, 
the dominant contribution to the difference is from 
the #24 term omitted in [30] which yields a positive 
contribution. At temperatures above 3000 °K, the 
dominant contribution is from the omitted xn and ^23 
terms which yield large negative contributions. In 
addition to these, the high value they used for x^ 
(larger by a factor of two) tends to reduce the differ- 
ences at temperatures below 3000 °K, but magnifies 
them above. The same general comments apply to the 
JANAF [31J comparison as seen in curve #1, except 
that the difference goes negative at somewhat higher 
temperatures and the positive deviations are consider- 
ably larger. Incidentally, the JANAF and Harrison and 
Kobe [28| (curve #4) calculations are nearly identical 
for temperatures up to 1000 °K. Curve #3 applies to 
a modified rigid rotator harmonic oscillator calculation 
as given in an older (now superseded) JANAF table. 
The comparison for curve #2 [29 1 illustrates the sen- 
sitivity to the omission of the hierarchy of states 
corresponding to large vi quantum numbers (vz > 3). 
The differences plotted in figure 2 exhibit essentially 
the same behavior as those in figure 1 except that they 
are somewhat larger. 
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FIGURE 2. Dimension/ess heat capacity differences. 



It is a pleasure to thank Joseph Hilsenrath for many 
helpful technical discussions, particularly relating 
to the use of OMNITAB to program the calculation. 
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